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Abstract 

Some aspects of a(j —fo mixing effects in the reaction pn — > da^ with perpen- 
dicular polarized proton beam are discussed. An angular-asymmetry param- 
eter A is denned to study those effect. It is shown that, for energies close to 
the production threshold, the angular-asymmetry parameter A(9, <p) is pro- 
portional to the a[] — fo mixing amplitude for arbitrary polar and azimuthal 
angles 9 and (p of the outgoing oq meson. This statement is also valid for 
arbitrary energies, but then only at polar angles 9 = 0° and 9 = 90°. The 
mass dependence of the differential cross section da/dm^o^ in the reaction 
pn — > dir°n in the presence of a[j — fo mixing is also discussed. 
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I. INTRODUCTION 



The nature of the lightest, virtually mass-degenerate, scalar mesons a (980) (I G J PC = 
1~0 ++ ) and /o (980) (0 + ++ ) is an important and still unsolved problem of hadron physics. 
The quark structure of these mesons is not well established at present time. This issue is 
also closely related to the very interesting problem of ao-/o mixing. A dynamical mechanism 
for this mixing close to KK threshold was suggested around 20 years ago in Ref. Q. Since 
that time a number of papers have been published in which different aspects of a$ -fo mixing 
and the possibilities to measure this effect have been discussed, see, for example, Refs. 0-0. 
Recently, in Ref. ||, some new arguments were presented in support of a fairly large a° — / 
-mixing intensity. Analyzing the experimental data [j9[|10| on the exclusive production of the 
ao(980), a2(1320) and /o(980) resonances in pp collisions at y/s = 29.1 GeV, the authors of 
Ref. H came to the conclusion that 80 ± 25 % of all a[j mesons come from the /o(980). 

The results of Ref. [§] provide certainly a fairly strong indication for a large a[j — / 
mixing. However, for more solid conclusions, in particular about the quantitative value of 
the mixing, a thorough analysis is needed. In this context it is important to keep in mind 
that the mixing can manifest itself in many different physical processes and observables. 
Thus, in this paper we want to discuss another consequence of the a[j —fo mixing effect, 
namely the angular asymmetry of the reaction 

pn-^da° /f . (1) 

with unpolarized and polarized proton beam. Furthermore we study the influence of the 
finite width of the scalar mesons on the unpolarized as well as the polarized differential cross 
section. All these investigations are performed with special emphasis on how to extract 
informations on the ao//o mixing. 

In the present paper we will address the do and /q mesons as resonances. However, it 
should be stressed that all the results derived in this paper do not need any assumption about 
the nature of those scalar mesons. Even if these scalar mesons are dynamically generated 
|TT| , [T2"[ |, in the proximity of the pole the scattering matrix can still be well approximated 
by the propagation of quasiparticles corresponding to elementary fields. The true nature of 
the propagating object is then hidden in the effective parameters of those elementary fields. 



This was shown rigorously by Weinberg for the case where inelastic channels are absent p3 

The paper is organized as follows. In Section 2 we review some common properties of the 
production of unstable particles. We start out from the case of a single particle and then, 
in Sect. 3, generalize the formalism to the two channel case relevant for the present study. 
Section 4 contains a discussion of the general structure of the primary production amplitude 
for small excess energies. Section 5 is devoted to the study of different observables. The 
main emphasis is put on the irr] final state. Two differential observables are analyzed in 
detail, namely the double differential cross section d?cr/dQ/dm 2 as well as the analyzing 
power. The former observable was studied recently in a series of papers, however, without 
investigating the m 2 dependence more thoroughly. The latter, on the other hand, was not 
discussed in detail before at all. We close with a short summary and conclusions. 
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II. THE PRODUCTION OF UNSTABLE PARTICLES 



Let us start with the general expression for the cross section of the reaction NN — > dX , 
where X denotes the decay products of an unstable meson. We start the discussion by 
assuming the presence of one mesonic resonance only. The generalization to more than one 
resonances will be done in the next section. The main purpose of this section is to introduce 
our notation and to remind the reader on the impact of a finite decay width on observables. 

The reaction cross section in the center-of-mass (CM) system be can expressed as 

da x = (2?r) 4 — - \A x \ 2 d$^ +1 (P; p d , k h . (2) 

4 lPlv s 

Here p denotes the CM momentum of the initial protons and P the total initial four- 
momentum. The reaction amplitude A x is given by 

A x = W x (k 1 ,...,kJG(m 2 )M . (3) 

where the primary production amplitude is denoted by M.. We assume the unstable meson, 
whose propagation is described by G(m 2 ), to decay into the /i particles of the final state 
X through the vertex function W x . m 2 is the total invariant mass of the final /i-particles 
system and is given by m 2 = (J^ki) 2 , with k iy i = 1, \i being the four- momenta of the \x 
decay particles. The phase space of the final fi particles and the final deuteron (with the 
four-momentum p d ) is defined as 

d^ +1 (P; Pd , h, fc M ) = S\P - Pd -k) J] t=^- 

= {2Ttfdm 2 d§ 2 {P- p d , k) dQ^k; h, fc M ) (4) 

where E d and uji are the energies of the final deuteron and the z-th decay particle, respectively, 
and k = J^ki- The latter recursive formula allows to study the propagation and the decay 
of the unstable particles independently of the production mechanism itself. 

It is convenient to introduce what we would like to call partial spectral functions p x , 
which is given by 

p x (m 2 ) := (2vr) 3 | d§^(k; k±, k^) \G(m 2 )W x \ 2 . (5) 
Note that unitarity demands that 

Y / Px(m 2 ) = --lmG(m 2 ) , / dm 2 ^ Px(m 2 ) = 1 , (6) 
x 71 J x 

where the sum runs over all decay channels of the unstable particle that are open at a certain 
given value of m 2 . 

Using Eqs. (|]) - (||) we can write 

*°* 1 /fi nA1 | W ), (7) 



dm 2 dVtk 647r 2 \\p\ I s 
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where k is the relative momentum of the meson system X with respect to the deuteron. 
Naturally, \k\ is a function of m 2 , namely 

where denotes the mass of the deuteron and the function A is given by 

X(x,y,z) = (x - y - z) 2 - Ayz. (9) 

In this paper we neglect a possible final state interaction (FSI) of the produced mesons 
with the deuteron (this issue was discussed for the reaction pp — > da^ in Ref. |T4]|). Then, 
close to threshold, the energy dependence of the matrix element M. is determined by the 
lowest possible orbital angular momentum I of the meson system X with respect to the 
deuteron. Specifically, it will be given by k l with k = \k\ being the corresponding relative 
momentum. 

Let us now consider the case that interests us in the present study, namely the decay of 
a scalar resonance or quasiparticle into two pseudoscalar mesons. For any given interaction 
of the scalar mesons with the pseudoscalars, the physical two point functions for the scalar 
mesons can then be constructed by solving the Dyson equation, where the self-energy S is 
given by the leading order two pseudoscalar loop diagrams. The equation is illustrated in 
Fig. |]. The strongest energy dependence of the loops is introduced by the unitarity cut 
and its analytical continuation below the threshold. The remaining piece can be assumed 
constant and can be absorbed in the physical mass Mr. Under these assumptions the 



physical propagators of the scalar mesons can be described by a Flatte form [15], namely by 



X 



where 



r 



X 



\WxV\kx 
8ttM| 



with kx being the CM momentum of the meson system X. Below the threshold of the 
production of the final state X the analytic continuation of Tx is to be used. 



III. THE MIXING AMPLITUDE 

We now turn to the case that is relevant for the present paper, namely when there are 
two propagating unstable particles with different isospins. Then the self-energy £ exhibits 
a matrix structure in isospin space. If we now allow the propagating particles to mix then 
the self-energy matrix develops non-diagonal elements. Note, since the mixing particles a 
and fo are essentially mass degenerate all the mixing should be completely dominated by 
its effect on the meson propagation itself. Along the lines of this reasoning we completely 
neglect any isospin violation in the primary production amplitude Ai and in the coupling 
of the pseudoscalar mesons to the resonances. In short, we assume the resonances to be 
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produced in pure isospin states. The physical states, however, are those that propagate. 
The physical states of a and f are consequently no longer pure isospin states. 
We start from the inverse propagator in the isospin basis, 

G^V 2 ) = 

where mi and denote the bare masses in the isospin 1 and channel, respectively. The 
diagonal elements of the self-energy are given by 

e u = vl\ k + (ii) 

EOO — ^°KK + ff 7T7T (12) 

and the isospin breaking is generated by 

Eio = E i = a + Hk+k- — N-k°k° • (13) 

In these formulas U xy denote loops integrals with particles x and y propagating. The loops 
are assumed to be renormalized - the renormalization constants are either absorbed into the 
physical masses of the scalar mesons or into the complex quantity a (see also the discussion 
of this question in Ref. [|16j)- The loops denoted by U.g K are to be read as containing 
contributions from charged kaons as well as neutral kaons, coupled to the isospin / as 
indicated in the superscript. It is also implied that the physical masses are used in those 
loop integrals. The mass difference between the charged and the neutral kaons is the source 
for the isospin breaking generated by the loops of Eq. ( p~3|) . Without this difference the two 
kaon loops would cancel each other. The complex quantity a contains the effect induced 
by tt-t] mixing p7|, as for instance discussed in Ref. M, as well as a possibly existing direct 



a o~fo transition. Naturally, the imaginary part of a, induced by the 7r°— 77 transition, can 
be easily estimated as was done in Ref. ||. 

The propagation of the physical states is described by the eigenvalues of G, i. e. 

G ao = 2/(g u + g 22 + r) and G fo = 2/(g u + g 22 ~ r) , (14) 



where r = J(gn — 922) 1 + 4g 2 2 . The physical masses of the scalar mesons a and f are 



given by the real part of the zeros of Eq. ( \A ) . 

As usual the full propagator G can be expressed in terms of the physical propagators G ao 
and Gf plus a parameter for the mixing, which can be identified as 5 = 2g 12 / (gu — g 22 ) . In 
the case of p-u mixing or tc-t] mixing 5 is guaranteed to be small due to a large denominator 
9n — 922- In the former case there is a large difference in the widths whereas in the latter 
the masses of the mixing particles are rather different. However, the experimental evidence 
that is available for the masses and widths of the scalar mesons ao and /o suggests that 
here both these quantities might be very similar ||18|| . Thus, even gu = g 22 is not excluded. 
Also, we want to emphazise that 5 can not be interpreted as a mixing angle because it is 
necessarily complex valued. 

For convenience we write the primary production amplitude M. as a vector in isospin 
space, 
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M 

and the final production vertex Wx as a matrix 

W - 

Note that here is the vertex for the coupling of the irr] system to an isospin-one particle 

as given by an elementary Lagrangian, whereas is the vertex where the irr] system 

mixes first into the tttt system which then couples to an isospin-zero particle. We thus find 
for the complete production amplitude of the final states X = (r/n) or (nn) (c.f. Eq. (Q) for 
the one channel situation) 

A = WtQM = Wt -( 922 ~ 9l2 )M = W^G ao G fQ ( 922 ~ 9u )m . 

gn922 - #12 V ~9n gn ) \ -912 911 J 

Basically all investigations on the a[j — f so far don't use the propagators of the physical 
states, i.e. they don't use the eigenvalues of G. Rather they assume that the ao and /o 
mesons still propagate as proper isospin states. This leads to expressions that are very 
easy to interpret and which are illustrated in Fig. 2 for the case of X = (irrj). Within the 
formalism developed above this assumption corresponds to the approximation 

Gf G ao = —— -j- — — — — = G fo G ao , (15) 

9\\922 — 9\2 9\\922 

where Gf and G ao are now propagators of these mesons in the isospin and 1 states, 
respectively. Indeed, the results presented in Ref. suggest that 

g 2 i2 <10 - 2 



911922 



i.e. that the approximation ([IB]) is not unreasonable. 

Since the main focus of the paper is the irr] channel it is convenient to introduce a mixing 
parameter £ by 

^(m) = S 01 (m)G' /0 (m 2 ), (16) 

where S 01 was defined in Eq. flTB| ) and is related to the a° — / transition amplitude. With 
this definition we get for the transition amplitude for the reaction pn — > dnrj 

= W^G^iM! + £(m)M ) + G fo M + wfkWG^Mi . (17) 

The last term involves mixing in Wt v) as well as in £(m) and is therefore a second order 
correction. Close inspection reveals that in the region m ~ 2m^ the third term in Eq. (|TT 
should be suppressed by a factor 



gaon^oiim^ - ml) 



0.12 
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compared to the second term. In this estimation we used the values Eoi ~ 5000 MeV 2 
determined from KK decay loops (see Fig. 2 in Ref. ||) at m ~ 2rriK and the irn mixing 
amplitude A^o^ ~ —5000 MeV 2 (see Ref. and references therein). We also use the masses 
and widths m ao = rrif Q = 980 MeV/c 2 and r ao = T/ = 50 MeV/c 2 and coupling constants 
g ao -Kri and gf 07in that are determined from the equations r ao = g 2 ^^ / '(87rm 2 ) and T/ = 
3(7 2 07r7r g 7r7r /(167rm 2 o ), where q nrj and q nn are corresponding relative momenta. Consequently, 
both terms - the third and the last in Eq. flT7p - will be neglected in what follows. 



IV. STRUCTURE OF THE PRIMARY PRODUCTION AMPLITUDE 



As we argued in the former section one can safely assume the primary production am- 
plitudes Aij as isospin conserving. In this section we demonstrate how to construct the 
effective interaction relevant for the transitions pn — > d + (scalar meson) in the close-to- 
threshold regime. 

As pointed out in Refs. |L9|||, if the scalar meson has isospin 1 it can only be produced 



in a p wave. Since the initial pn system has to be in an isotriplet state and has to have odd 
parity, the Pauli principle requires that it is also in a spin-triplet state. Thus, the effective 
transition operator for the isovector final state has to be linear in k and has to have an odd 
power of p (hereafter k and p denote the final and initial relative three-momenta in the 
considered reaction, respectively). It also has to be linear in both the spin S := <f)J<j2c4>2 
of the initial nucleons pair and the polarization vector e of the outgoing deuteron. These 
constraints lead to the reaction amplitude of the following type 

Mi = a(p-S)(k- e*) + b{p-k){S-e*) + c{k- S) (p -e*) + d(p- S) (p • e*) (k ■ p) . (18) 

The coefficients a, b, c and d are independent scalar amplitudes. They may depend on the 
total CM energy and are necessarily complex since they contain the initial state interaction. 
For the isospin 1 amplitude relevant here the phase induced by the two nucleon unitarity 
cuts can be related to the NN phase shifts, as pointed out in Ref. pOfl . However, for 
the isoscalar initial state (required for the production of the f meson, cf. below) there is 
no phase-shift analysis available at the relevant energies due to a lack of high energy pn 
scattering experiments [pi ]. Therefore, in this paper we will not consider the effects induced 



by the initial state interaction. Let us mention, however, that their effects do not influence 
the qualitative aspects discussed here. Note, that the scalar amplitudes a, b, c and d are 
expected to have quite smooth energy dependence and, accordingly, can be considered as 
constant in the near threshold region. 

In Ref. |14]| the transition amplitude was given in a partial wave decomposed form. Ob- 



viously the two descriptions are equivalent, however, we regard Eq. ( |I8"D as more convenient 
for the construction of polarization observables. 

The amplitude relevant for the near-threshold production of an isoscalar scalar meson 
can be constructed along the same lines and is given by 

M = f(S-e*)+g(p-S)(p-e*). (19) 

Similar expressions for the amplitudes (|18"D and (|19|) can be found in Ref. [f7p2|. Note, 
however, that the terms proportional to g and d were omitted in that work. 
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It is obvious from Eqs. fllSl) and (|I~9"D that .Mi oc k and A^o oc const, near the threshold. 
For the isospin conserving case we therefore get 



; (pn — > cZa ) oc Q 3 ^ 2 and -j-^j(pn — > d/o) oc Q 1 ^ 2 



c?m 2 " <im 2 

where Q(m) = y/s — — m is the excess energy. Thus, if we study the production of an 
isospin 1 final state, namely n°n, the mixing with the isoscalar (/ ) state will be kinematically 
enhanced [HI. 



V. STUDY OF OBSERVABLES 

A. Invariant mass spectrum 

Naturally the simplest observable is just the invariant mass spectrum of the reaction 

pn — > dX. As long as we look at energies/values of m 2 such that Q(m 2 ) is small, the 

spectrum should show a resonant peak ~ J^Px reflecting the m 2 dependence of the fo 

x 

propagator, cf. Eq. @. It is possible to deduce the branching ratios of a particular 
resonance by a fit of the Flatte distribution to the mass spectrum. This was demonstrated, 



e.g., for the case of the czq in Ref. [23]. Since so far little is known about the fo this very 
easy to measure observable is very interesting. In addition, the number of all events in this 
peak is proportional to 

\Mo\ 2 = 3\f\ 2 + 2Re(fg*)p 2 + \g\ 2 p\ 
cf. Eqs. (U), (0) and ([19]). Note that the knowledge of the amplitudes / and g is important, 



if we want to deduce quantitative informations about the /o/ a o mixing amplitude from the 
reaction pn — > dX. 

Definitely the most interesting observables are those concerning the channel pn — > dnn, 
since here the isospin conserving amplitude enters in a p wave whereas the isoscalar ampli- 
tude, that can be mixed in via isospin violation, enters in an s wave. This is the kinematical 
enhancement mentioned at the end of the last section. 

Let us discuss the mass distribution da/dm^^ for the 7r°n system in the reaction pn — > 
dir°n. Note that near the threshold of the reaction pn — > da® this mass distribution should 
be sensitive to the magnitude of a[j — fo mixing as well as to the mixing mechanism. 

Consider first the limiting case of isospin conservation, so that the a(j — fo mixing is 
absent. In this case the mass spectrum of the ir°n system should coincide with the spectrum 
of the ir + r] system from the reaction pp — > rfa^ - apart from an overall (isospin) factor of 
0.5 0[. The amplitude of the reaction pn —>■ dir ri is then given by the diagram of Fig. 2a. 
In this case the invariant 7r°?7-mass distribution is determined by the phase space and P( 7rr) ), 
which was defined in Eq. (§), and is dominated by the a propagator (cf. Eq. (|10"D). The 
7r°ry-mass distribution obtained under these asumptions, is shown in Fig. 3 by the dotted 
line. The calculations are done at T p = 2645 MeV/c where experimental results can be 
expected from the ANKE collaboration soon [^] . The parameters used here are: M ao = 980 
MeV/c 2 , r ao = 50 MeV/c 2 . The mass-dependent width is described by a Flatte-like form 
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r ao (m) = r ao + T K ^(m), which takes into account the KK decay channel (see Eq. (8) in 
Ref. H). As we discussed in the previous Sections, the isospin conserving production of the 
a[J meson can take place only in a p wave and hence is suppressed near threshold. 

If isospin is not conserved and afj — fo mixing is present, then the ttt] system can be 
produced also in an s wave through the fo- The corresponding diagram is shown in Fig. 2b. 
For this mechanism the 7r r] mass spectrum is determined not only by the a[] propagator, 
the partial wave of production and the phase space but also by the fo propagator and by 
the mass dependence of the nondiagonal self-energy matrix element S i(m) (cf. Eqs. ( |13|) 
and QTD): 

£(m) = £ i(m)G /o (m 2 )~ 



m 2 -M 2 o + iM fo T fo (m 2 



For illustrative purposes we will use M/ = 980 MeV/c 2 and the same Flatte-like width that 
was used for the ao meson above, i.e. Tf (m 2 ) = r ao (m 2 ). 

Of course, now the invariant mass plot becomes sensitive to the specific mixing mecha- 
nism. For example, if we assume that the leading contribution to Soi(m) comes from the 
direct a[] — fo transition mechanism, i. e. the term a in Eq. (|T3"D, then there should be no 
drastic dependence on the invariant mass m of the 7r°rj system. The 7r°?7-mass spectrum, 
corresponding to the mechanism of Fig. 2b and with Eoi(m) = const, is shown in Fig. 3 by 
the dashed line. (All the distributions are normalized to 1 at the maximal values!) On the 
other hand, if S i (m) is determined mainly by the KK loop (the term U K + K - — IT^-o^o in 
Eq. ([IB])), it should be strongly enhanced near the KK threshold (see, e.g., Refs. |],|4|.|(| ) . 
The corresponding 7r°?7-mass distribution is shown in Fig. 3 by the solid line. 

Thus, we conclude that the n°r) mass spectrum is rather sensitive to the a® — fo mixing 
mechanism. That is why the study of da / 'dm^o^ in the reaction pn — > dir n should allow to 
shed light on the nature of the lightest scalar mesons ao(980) and /o(980). 



B. Analyzing power for the process pn — > dirij 



1. Near threshold 



In our approximation the total amplitude Ai of the resonant 7rr/-pro duct ion with a[j —fo 
mixing effects taken into account may be written as (cf. Eq. (|17|)) 



M = Mt+^Mo, 



(20) 



where the amplitudes Aii and M.o ar e given by the Eqs. fll8|) and (|T9|) and the mixing 
parameter £ was defined in Eq. (IB). 



Let us introduce the polarization of one of the nucleons C, through £ 



>I Or (prfj 



(1 + £-<t)/2. The matrix element (|20|) squared and averaged (summed) over the polarizations 
of the initial neutron (final deuteron) is then given by 



\M\ 



(H 2 + |c| 2 ) p 2 k 2 + 3\B\ 2 + \D\Y] + {p-k) ReH 
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+ p 2 Re [a*D (p ■ k) + B*D] + (£ • [k x p])ImiZ\ 



(21) 



where 

B = b(p-k) + f f, D = d(p-k)+g£ ) H = a*B + a*c(p ■ k) + B*c + D*cp 2 . 

Let the polarization vector £ be directed along the x-axis (£ J_ p). Then, in terms of the 
angles 9 and (p, we have p ■ k = pk cos 9 and £ • [fe x p] — (pk sin 9 simp (( = Using 
those expressions Eq. (|2lf) may be written as 

(22) 



where 



\M(9,(p)\ 2 = C + d cos 9 + C 2 cos 2 9 + ( sin9 sin p (D + D 1 cos I 



Co = l(\a\ 2 + \c\ 2 )p 2 k 2 +(\f\ 2 + l -\f + p 2 g\ 2 ) |e| 



Ci =pfcRe([(a + 6 + c + p 2 rf)*(/ + p 2 (7) + 267] £ 



2 ;„2 



C 2 = 



|6| s 



|6 + p 2 d| 2 + Re (a*c + (a + c)*(6 + p 2 d) 



(23) 



J D =pA;Im([a7-c*(/ + p 2 (?)J £j , 
D 1 = p 2 k 2 lm(a*b + a*c + b*c + p 2 d*c) . 

We can now generalize the definition of the asymmetry A given in Ref. to the polarized 
situation. Using the short hand notation 

d 2 a 



a(m; 9, (f) = 



dm 2 dVt 

we now define the angular-asymmetry parameter A through 

cr(m; 9, ip) — cr(m; n — 9, ip + 7r) 



(24) 



cr(m; 6 1 , <p) + cx(m; 7r — 6 1 , + it) 
Ci cos 6> + ( D sin sin <^ 



Co + C 2 cos 2 9 + ( Di cos 9 sin 6* sin </p 
It follows from Eqs. (||) and (g|) that Ci = D = for £ = and thus 

A(m; 9, (p)^ =0 = 



(25) 



(26) 



when isospin is conserved. 

Let us discuss some specific features of the expression (p5"D in the following. First note 
that there are two different terms. The first term does not depend on the polarization ( of 
the proton beam and on the angle ip. Therefore, Eq. (|25| ) implies that 



A( m -9 = 0°,<p) 



Co + C2 



(27) 
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This particular result for the asymmetry A was derived earlier in Ref. || but not in terms 
of the amplitudes a, b, c, d, f and g. The discussion in Refs. [0,0], on the other hand, takes 
into account only the amplitudes a, b, c, and /, but not the amplitudes d and g. 
At 8 = 90° and for ( ^ we get 

A( a Qn o n , D ° ■ 2(pklm([a*f~c*(f + p 2 g)] flsiny 

A(m;8 = 90 , <p = C 77- siny? = . . 28 

C (|a| 2 + |c| 2 )p 2 fc 2 + (2|/| 2 + |/ + p 2 g\ 2 ) K| 2 

Obviously the isospin-breaking effect for different angles 8 and <p depends on different 
combinations of the basic amplitudes. Thus, experimental information on the asymmetry 
obtained with polarized beam could allow to deduce additional constraints on the amplitudes 
a, b, c, d, f, and g of the a - and f production reactions. With regard to that let us mention 
that the model of ao-meson production discussed in Ref. || , which was based on the impulse 
approximation, leads to A(m; 8 = 90°, tp) = - because in that model a = c = 0. 

Throughout this paper we make the assumption that isospin breaking takes place only 
in the propagation of the scalar mesons. This implies that the mixing strength is the same 
for all partial waves. On the other hand, if a significant mixing takes place also in the initial 
NN interaction or in the primary production amplitude A4i, there is no reason to expect 
this mixing as being partial- wave independent. Thus the study of polarization observables 
allows to examine this basic model assumption in a clean way. 

In addition it will be very interesting to study the m dependence of A(m; 8, <p) in detail. 
By construction A projects on the isospin breaking pieces of the amplitude. Therefore, the 
m dependence of A gives direct access to the m dependence of S i(m) and thus to the 
mixing mechanism. 



2. Higher energies 

As mentioned above the asymmetry A(m; 8, if) (Eq. (|25|)) will become zero if a[j — /o 
mixing is not included, i.e. for £ = 0. This result is in line with the observation made 
in Ref. 0] that the asymmetry A(8) is identical to zero in the reaction pn —> dn for 



the unpolarized case. Let us emphasize, however, that the vanishing of the asymmetry 
A(m; 8, ip) in the limit £ = for the polarized case is true only near threshold, where the 
ao meson is produced in a p wave and the fo meson in an s wave, see Eqs. fll8|) and (fL9|) . 
The expressions (|T8| ) and ( [T9D do not take into account contributions from higher partial 
waves, specifically they do not include, e.g., the d wave production of the ao meson and the 
p wave production of the fo meson. In the following we want to consider this question in 
more detail, assuming that isospin is conserved. 

The amplitudes for a - and fo production in the reaction ([TJ) can be written in the most 
general form as 

Mj = cf>Ja 2 (Fj + Gj-a)^, (29) 

where <pi t 2 are the spinors of the nucleons. The two terms Fj and Gj correspond to SWat — 
and SVat = 1 , respectively, where Snn is the total spin of the initial iVW system. The 
subscript I denotes the isospin of the produced scalar meson. The scalar and vector functions 
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Fj and Gi depend on the vectors p, k and e. Let Lnn and Lud be the angular momenta 
in the initial NN and in the final meson+(i systems, respectively. It follows from the 
conservation of quantum numbers together with the required antisymmetry of the system 
with respect to the initial nucleons that the F\ {G\) term in Eq. ( p9|) should contain only 
the contributions from even (odd) values of Law and Lud- One can also see that the F 
(Go) term should only contain the contributions from odd (even) values of Law and LMd- In 
the isovector (I = 1) case the s wave (Law = L>ud = 0) is forbidden, i.e. the contributions 
to the Fi term start with a d wave. For the case where one of the nucleons is polarized, the 
expression for the squared matrix element (^) is given by the form 

1 



l-M/l 2 = \ [l^/l 2 + ({FiG} + F/Gj] • C) + (G*i ■ G T ) + i (C • [Gj x G}})\ , (30) 



where £ is the proton polarization vector. The second term in the expression fl30|) corresponds 
to the interference between the Lj and Gj amplitudes. Thus, this term is an odd function 
of the final momentum k. All other terms in Eq. (|30| ) are even functions of k. We see that 
the second term is responsible for the angular asymmetry in the reaction ([]]) even if isospin 
is conserved. It vanishes in the case of unpolarized protons, i.e. A(m; 9, <p) = for C = 0. 

Note that the amplitudes (|T8|) and ( |T9"D are special cases of the general form Eq. (|29"D with 
Fj = 0. Also, setting Fj = in Eq. flBTf ) we always get A(m; 9, (p) = if isospin is conserved. 
The amplitudes ( fL8|) and (|T9|) correspond to the Gi terms in Eq. ( |30| ) in the lowest-order 
approximation with respect to the final momentum k, i.e. keeping only contributions that 
are at most linear in k. The next-order term with respect to k in the amplitude for ao-meson 
production is the <i-wave contribution. This term (of the order ~ k 2 ) contributes to F\ in 
Eq. (p9|) and modifies the amplitude M.\ given in Eq. (fOp, i.e. 



Mi -> A^i + Fx 0fa 2 2 , F 1 =e(e*- [p x k]){p ■ k) (31) 

where e is a scalar amplitude. It is clear that the asymmetry A(m; 8, ip) (|25|) , calculated with 
the amplitude (|3~ij), will get non-zero contributions ~ k 3 for small k due to the interference 
between p and d waves. It follows from Eq. ([31]) that F± = at sin 6* = or cos = . 
In general, the Fx term contains all contributions with even values L NN ,L Md = 2,4,6, ... . 
Thus, each contribution should include the factor e* ■ [p x k] as well as the factor p ■ k. 
Therefore, even in the general case Fx — if sin 6* = or cos# = 0. Accordingly, for the 
case of conserved isospin we get 

A(m; 9 = 0, cp) = (32) 

at any k - and not only near threshold (compare this result with that from Eq. (|26|)). This 
result p!3) confirms a general statement made in Ref. on the absence of a forward- 



backward asymmetry for the reaction pn — > dX° in the limit of isospin conservation. 

It is remarkable that in addition to the condition (^) for the case of polarized proton 
beam we get also the nontrivial result that 

A(m; 9 = 90°, <p) = (33) 

if isospin is conserved. Therefore, any deviation of A(m; 9 = 90°, ip) from zero is a direct 
indication for a° — fo mixing. Note that this statement does not depend on the number of 
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partial waves taken into account, i.e. it is valid for any excess energy and not only near 
threshold. Consequently, a measurement of the asymmetry at 6 — 90° should provide us 
with evidence on the a[j — fo mixing amplitude. We want to remark that this information is 
to be considered as complementary to the one that can be obtained with unpolarized beam. 
It should be mentioned that the vanishing of the asymmetry at 9 = 90° for the reaction 
pn — > da,Q does not follow from the theorem formulated in Ref. [24| and is new. 

Finally, let us discuss the situation at 9 ^ 0° and 9 ^ 90°. If the amplitude for a[j- 
meson production is taken of the form (|20|), i.e. isospin-breaking effects are included, then 



the leading-order contribution to the asymmetry A(m; 9, ip) are of the order ~ £/c (see 
Eqs. ( |23"D and (^5|)). The next-to-leading terms in the amplitude .Mi (see Eq. (^TJ)) give rise 
to contributions of the order ~ k 3 to the asymmetry from isospin-conserving terms. Thus, 
when studying isospin-breaking effects over a wide region of 9, one should consider a^-meson 
production in a rather narrow region of relative momenta k, i.e. at small excess energies, 
in order to suppress contributions to A(m; 9, ip) from isospin-conserving amplitudes with 
higher angular momenta. 

VI. SUMMARY 

Let us summarize the main results of this paper. The most sensitive observables for 
examining the a° —fo mixing amplitude are: 



i) the angular asymmetry A(m; 9, ip), as defined by Eq. (|25|) of this paper 



ii) the distribution of the effective mass of the it n system near the a° threshold 

With regard to the angular asymmetry A(m; 9, (p) for the channel pn — > d7r°n, it is expected 
to be rather large, i.e. in the order of ~ 10 _1 . For the case of an unpolarized beam the 
forward-backward asymmetry was estimated earlier in Ref. ||. The study of the asymmetry 
A(m; 9, ip) with perpendicular polarized proton beam should shed additional light on the 
mechanism of a>o production in the reaction pn — > da®. As we have shown the various basic 
amplitudes for this reaction give different contributions to A(m; 9, ip) at different angles. 

Since by constuction A projects on the isospin breaking pieces of the amplitude, the m 
dependence of A gives direct access to the m dependence of S i(m) and thus to the mixing 
mechanism. 

Note, in experiments with polarized beams isospin-breaking effects give contributions 
~ £&; to the asymmetry A(m; 9, ip) in leading-order with respect to the final momentum 
k. Isospin-conserving amplitudes give contributions ~ k 3 due to d-wave terms. Thus when 
studying isospin-breaking effects over a wider range of 9 and <p, one should consider a[]-meson 
production at small momenta k, or small excess energies. Accordingly, when studying the 
reaction pn — > d7T°n near the threshold of a® production, one should limit the invariant mass 
m of the 7i°r] system in a narrow region near threshold. Note that the contribution of the 
c?-wave to the reaction pn — *> dao may be monitored by a simultaneous measurement of the 
differential cross section. 

If isospin is conserved, then asymmetry effects should be absent at 9 = 0° or 9 = 90°, 
even in the polarized case. The angle 9 = 0° is not conclusive for the reaction with polarized 
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beam, since all polarization effects are proportional to £ sin # cosy? anc l vanish at sin# = 0. 
Thus, the asymmetry for 9 = 0° coincides with the one for unpolarized nucleons. The case 
6 = 90° is much more interesting when studying the isospin-breaking effects in the reaction 
with polarized protons, since polarization effects are maximal at this angle. Anyway, in 
either case (6 = and 9 = 90°) a non-vanishing asymmetry A(m; 8, (p) can come only from 
isospin-breaking effects. 

Very important informations on a° — fo mixing can be also expected from a study of 
the mass distribution da / dm n o v for the 7r°?7 system in the reaction pn — > dir !]. If the 
isospin is conserved then the mass spectrum of the 7T°rj system from this reaction should 
coincide with the one for the n + 7] system from the reaction pp —>■ da^ . Differences in the 
mass distributions of the n r] and the 7r + ?7 systems are a clear signal that the isospin is not 
conserved and that a[j —fo mixing is present. The mass distribution da/dm^o^ for the 7r°?7 
system is also extremely sensitive to the mechanism of a[j — fo mixing. For example, in the 
case where a§ — fo mixing is dominated by the KK loop, the n°r] mass spectrum will be 
strongly enhanced near the KK threshold. 

Thus, we conclude that the study of the reaction pn —>■ da^ is extremely useful to under- 
stand the mechanism and the magnitude of the a{] — fo mixing amplitude. The knowledge 
of this mixing amplitude should allow us to shed light on the nature and the quark content 
of the lightest scalar mesons a (980) and / (980). 

Several interesting observables where not discussed in this paper and are to be studied 
in a subsequent work. One example is the differential cross section for pn — > dK + K~ . Since 
this final state does not have definite isospin the differential cross section develops a forward- 
backward asymmetry at larger excess energies Q. Models predict that the kaon loops are 
an important source for the mixing Q. Therefore, the measurement of a forward-backward 
asymmetry for the process pn — > dKK (the sum of pn — > dK + K~ and pn — > dK°K°) 
should help to pin down the —fo mixing amplitude. Corresponding measurements of 
these decay channels should be feasible at the ANKE and TOF facilities, respectively, at 
the COSY accelerator in Jiilich. 

Finally let us consider the reaction dd — > 4 He7r°?7. Evidently, here isospin is not 
conserved. Near the ao threshold this reaction should proceed only through the chain 
dd — > 4 He fo —> 4 He —>■ 4 He Tr°n. Thus, the 7T°rj mass spectrum of this reaction should be 
even more sensitive to the ad — fo mixing mechanism, than the one for the reaction pn — > dad 
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FIGURES 



FIG. 1. Graphical representation of the Dyson equation. The thin solid line denotes the bare 
propagator whereas the thick solid line stands for the dressed/physical propagator. The sum 
appearing in the expression for the self energy is assumed to run over all decay channels X. 

FIG. 2. Different contributions to ttt] production in leading order in the mixing. Here /o — ao 
mixing is denoted by a circle and ir — r\ mixing is denoted by a cross. Note that /o — ao mixing 
can occur via KK loops as well as via direct mixing, as explained in Sect. 3. The numbers in the 
primary production vertices (denoted by a big circle) indicate the isospin of the relevant amplitude. 

FIG. 3. Spectra of the 7T rj invariant mass for the reaction pn — > dir°r] at T p = 2645 MeV/c. 
The dotted line corresponds to a reaction mechanism where an ao meson is produced in a p-wave 
and then decays into the TT°r] system. The dashed line corresponds to the production of an /o 
meson in an s-wave that mixes into an ao with Eio = const, cf. Eq. (|l^). The solid line shows 
results where Sio is evaluated from the KK mixing mechanism. Note that all the distributions 
are normalized to 1 at their maximal value. 
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